Temporal variations of solar rotation rate at high latitudes 
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ABSTRACT 

Frequency splitting coefficients from Global Oscillation Network Group 
(GONG) and Michelson Doppler Imager (MDI) observations covering the pe- 
riod 1995-2001 are used to study temporal variations in the solar rotation rate 
at high latitudes. The torsional oscillation pattern in the Sun is known to pene- 
trate to a depth of about O.li? with alternate bands of faster and slower rotating 
plasma. At lower latitudes the bands move towards equator with time. At higher 
latitudes, however, the bands appear to move towards the poles. This is similar 
to the observed pole-ward movement of large scale magnetic fields at high lati- 
tudes. This also supports theoretical results of pole-ward moving bands at high 
latitudes in some mean field dynamo models. The polar rotation rate is found to 
decrease between 1995 and 1999 after which it has started increasing. 

Subject headings: Sun: oscillations; Sun: rotation; Sun: interior 



1. INTRODUCTION 

The rotation rate of the solar interior can be determined by inverting observed rotational 
splittings of solar oscillation frequencies (Thompson et al. 1996; Schou et al. 1998). With the 
accumulation of GONG and MDI data over the last six years, it is now possible to study the 
temporal variation of the rotation rate in the solar interior. The rotation rate is known to 
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show temporal variations, with bands of faster and slower rotating regions moving towards 
the equator with time (Schou 1999; Howe et al. 2000; Antia & Basu 2000) similar to torsional 
oscillations observed at the solar surface (Howard & LaBonte 1980; LaBonte & Howard 1982; 
Snodgrass 1992). This pattern is found to penetrate to a depth of about 0.1f? Q . Torsional 
oscillations are believed to arise from nonlinear interactions between magnetic field and 
differential rotation. As such, they should provide a constraint on theories of solar dynamo. 
Covas et al. (2000) considered an axisymmetric mean field dynamo model to study temporal 
variations in rotation rate and magnetic field in solar interior. They find temporal variation 
in the rotation rate which shows pattern similar to torsional oscillations at low latitudes, with 
bands of faster and slower rotating regions moving towards the equator with time. But at 
high latitudes they find that these bands migrate pole-wards. As far as surface observations 
go, some magnetic features are seen migrating pole-wards at high latitude (Leroy & Noens, 
1983; Makarov & Sivaraman 1989). This pole- ward movement may be crucial for magnetic 
field reversal during the solar cycle. Thus it is of interest to check if the observed zonal flow 
pattern also moves pole-wards at high latitudes. Inverters thus far have largely ignored the 
variation of the rotation rate near the polar regions, where the measurement of the rotation 
rate is likely to be less reliable. Preliminary helioseismic work by Howe et al. (2001) and 
Basu & Antia (2001) suggests that that there could be a pole- ward flow at high latitudes. In 
this work, we make an attempt to study the temporal variations of the rotation rate near the 
poles. We first test the ability of inversion techniques to resolve the rotation rate near the 
poles using artificial data. We then apply these techniques to study the temporal variations 
in the rotation rate using observed splitting coefficients. 



2. DATA AND TECHNIQUE USED 

We have used data sets from GONG and MDI for this investigation. These sets consist 
of the mean frequency and the splitting coefficients for each (n,£) multiplet. We use the 
GONG data for months 1-55, which cover the period from 1995 May 7 to 2000 October 6. 
We use all available 54 data sets each covering a period of 108 days with a spacing of 36 
days between consecutive data sets. The first data set covers only 36 days (a GONG month 
being 36 days). The MDI data (Schou 1999) consist of 23 non-overlapping data sets each 
covering a period of 72 days, starting from 1996 May 1 and ending on 2001 April 4, with 
some gap in between when the SOHO satellite was out of contact. 

We use a 2-dimensional Regularized Least Squares (RLS) inversion technique to infer 
the rotation rate in the solar interior from each of the available data sets. The details of 
inversion technique are described by Antia et al. (1998). Since in this work we are mainly 
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interested in the rotation rate at high latitudes where the inversion technique are somewhat 
uncertain, we first perform a series of tests using artificial data sets to ascertain the reliability 
of the inversion procedure. For this purpose, we use a model rotation profile and calculate 
the splitting coefficients for the assumed rotation rate. Then we add random errors with 
same standard deviation as the estimated errors in observed data sets and perform inversion 
using the same regularization as that used for the real data. We construct artificial data sets 
using test profiles of the form 

n(r, 9) = 460 - 55 cos 2 9 - 55 cos 4 9 + A ex p(-(^^) 2 ) exp(-( C ° s2 - - C ° ) 2 ) (1) 

w r w c 

where Q(r,9) is the rotation rate in nHz and A, r , w r , Co, w c are suitably chosen constants. 
In Fig. 1, the inverted profiles are compared with the actual ones for a few of these sets using 
errors in MDI data sets. It is clear that inversions are able to reproduce relatively sharp 
peaks at high latitudes too, although the exact shape may not match the actual peak. The 
radial extent of the peak is reasonably well reproduced by the inversions. Thus the inversion 
results appear to be reliable even at high latitudes. 



3. RESULTS 

To identify the time varying component of rotation rate, we take the time average of all 
the results obtained from the different GONG (or MDI) data sets, and then subtract this 
mean from the rotation rate at each epoch to get the residual. This residual, SQ, contains 
the time- varying part of the rotation rate. Fig 2 shows the contours of constant residual as a 
function of time and latitude at a depth of O.O2i? below the solar surface as obtained from 
both GONG and MDI data. The MDI results have a gap during 1998-99 when no data were 
available due to problems with the SOHO spacecraft. In order to facilitate comparison of 
the inversion results with surface observations we show the rotation velocity = rQcos9, 
where 9 is the latitude. Also shown in the figure are the contours obtained for surface 
rotation velocity using Doppler measurements from Mt. Wilson (Ulrich 2001). For the figure 
with Doppler results, we show the north-south symmetric part of the rotation velocity since 
that is what the inversion results determine. In addition, we bin the Doppler data over the 
time intervals covered by the GONG data for a better comparison. In all these results we 
can see that bands of faster and slower rotation move towards equator at low latitudes. But 
beyond about 50°, the bands appear to move towards the poles. Theoretical results of Covas 
et al. (2000) based on a mean field dynamo model also show this feature. The latitude at 
which the transition from equator-ward to pole-ward movement takes place is also similar in 
their models. The pole-ward drift is clear in the GONG results which do not have any gap. 
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The surface rotation rate as measured using Doppler techniques (Ulrich 2001) also shows 
similar features, though the pole-ward migration of bands at high latitude is not as clear 
as in the helioseismic results, probably due to difficulties of making Doppler measurements 
at high latitudes. Helioseismic data represents an average rotation rate over the period of 
observations and this averaging smoothes out variations on short time scales in rotation rate, 
which have been observed in both Doppler measurements at solar surface (Ulrich et al. 1988; 
Hathaway et al. 1996) and in results from local helioseismic techniques (Basu & Antia 2000). 
To test the robustness of the pattern we have tried many different regularization parameters 
for inversion and all results show these bands. The GONG and MDI pattern will not match 
exactly as the temporal mean which is subtracted to calculate the residuals, is taken over 
different time intervals. 

In addition to having these bands, the rotation rate in the polar region has been de- 
creasing with time during the period 1995-1999. To show this variation clearly we show 
in Fig. 3 the rotation rate residuals at constant latitude as a function of time at a depth 
of O.O2i? below the solar surface. It is clear that there is a good agreement between the 
GONG and MDI results. At the latitude of 85°, which is the highest latitude for which we 
have tried to calculate inversion results, both GONG and MDI results show a clear decrease 
in the rotation rate residuals from 1995 to 1999, after which the residuals start increasing. 
The minimum in polar rotation rate is found to be in early 1999 in GONG data and slightly 
later in MDI data. At a latitude of 75° the pattern of variation is similar, but the amplitude 
of variation is much lower. The amplitude decreases rapidly as we move away from pole. 
It may be noted that in this figure we have plotted, the rotation rate, Q rather than the 
rotation velocity v^. The amplitude of v<j, variation does not increase as we approach the 
pole. 

To study the depth dependence of the changes in the polar rotation rate we show, in 
Fig. 4, the rotation rate residuals at a latitude of 85° at different depths. At r = O.98i? 
the errors are small and both GONG and MDI data show a clear time-variation with a 
minimum in early 1999 as discussed earlier. This panel also shows the North-South symmet- 
ric component of surface rotation rate as inferred from Mt. Wilson Doppler measurements 
(Ulrich 2001). It is clear that although the Doppler measurements have larger fluctuations, 
the basic trend is similar to the helioseismic results and in particular, these measurements 
also show the minimum at around the same period. The results for deeper layers have larger 
errors, but even so a similar variation is seen at r = 0.95R Q . At even deeper layers there is 
no clear temporal variation. Thus it seems possible that the depth to which the temporal 
variation extends is similar to 0.1i? Q found for the zonal flow pattern at low latitudes (Howe 
et al. 2000; Antia & Basu 2000). 
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4. CONCLUSIONS 

The solar rotation rate obtained from inversions of different sets of GONG and MDI 
data are used to study the time evolution of the rotation rate. The rotation-rate residuals, 
obtained by subtracting the time-averaged rotation rate from that at each epoch, show 
the well known pattern of temporal variation similar to the torsional oscillations observed 
at the surface, with bands of faster and slower rotation moving towards the equator with 
time at low latitudes. At high latitudes it appears that the bands move towards the pole 
instead, the transition between equator-ward and pole-ward movement appears to be around 
a latitude of 50°. Observations of magnetic features also show such pole-ward movement at 
high latitudes (see Leroy & Noens 1983; Makarov & Sivaraman 1989; Erofeev & Erofeeva 
2000; Benevolenskaya et al. 2001). Theoretical results of Covas et al. (2000) using a mean- 
field dynamo model also show this feature. Our inversion results therefore re-inforce the link 
between zonal flows and the solar magnetic cycle. Earlier works too have shown connections 
between zonal flows and the solar activity cycle (Antia & Basu 2000). 

The rotation rate in the outer layers of the polar regions varies with time, reaching a 
minimum in early 1999 after which it has started increasing. The time of minimum rotation 
rate at poles is distinctly before the maximum in solar activity. Similarly, it appears that 
the maximum rotation rate was achieved before the minimum in activity. These changes 
appear to persist till a depth of about 0.1f? Q , similar to the depth of penetration of zonal 
flow pattern at low latitudes. The observed temporal variations in the rotation rate should 
provide constraints on dynamo models. 
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Fig. 1. — Artificial data inversion results for several different model rotation profiles. In 
each panel the dotted lines show the true rotation rate, while other lines show the inverted 
profiles at different latitudes. The continuous line is for 85°, the short dashed line for 75°, 
the long dashed line for 67.5° and dot-dashed line for 60° latitudes. All the artificial profiles 
are for r = 0.95i? Q , w r = 0.02R Q , A = 20 nHz in Eq. (1). The panel (a) is for c = 1, 
w c = 0.1, panel (b) is for Co = 0.9, w c = 0.1, panel (c) is for c = 0.8, w c = 0.1, while panel 
(d) is for c = 0.8, w c = 0.05. 
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Fig. 2. — Contour diagrams of constant rotation velocity residuals at r = O.98i? obtained 
using 2D RLS inversion of GONG and MDI data. Also shown are the residuals for surface 
Doppler observation of Ulrich (2001). The continuous contours are for positive dv^, while 
dotted contours denote negative values. The contours are drawn at interval of 1 m/s. 
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Fig. 3. — The rotation-rate residuals at r = 0.98i? Q plotted as a function of time for different 
latitudes. The latitudes are marked in the figure. The results were obtained using 2d RLS 
inversion of GONG (left panel) and MDI (right panel) data. 
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Fig. 4. — The rotation-rate residuals at the latitude of 85° plotted as a function of time at a 
few selected radii as marked in each panel. The results were obtained using 2d RLS inversion 
of GONG and MDI data. The continuous lines show the GONG results with dotted lines 
showing the la error estimates, while the filled squares with error-bars show the MDI results. 
The thick dashed line in the topmost panel shows the results from Doppler measurement at 
the solar surface. 



